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1.  Introduction 


In  order  to  assess  the  characteristics  of  laser  propagation  above 

2 

the  sea,  it  is  necessary  to  estimate  the  distribution  of  C„  ,  the  structure 

N 

constant  for  index  of  refractions. 

The  index  of  refraction  depends  primarily  on  temperature  and  moisture; 
in  particular,  we  will  see  in  Section  3  that 
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where  C_  is  the  structure  constant  for  temperature,  B  is  the  Bowen  ratio 
and  r  is  the  correlation  between  high-frequency  temperature  and  moisture 
2 

fluctuations.  c  is  a  known  constant.  If  r  =  1,  this  equation  reduces  to 
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On  shipboard  only  relatively  simple  meteorological  observations  are 
usually  available,  e.g.  wind,  temperature  and  moisture  at  some  height  of 
order  10  m,  and  surface  temperature. 

2  2 

Therefore,  methods  have  been  proposed  to  estimate  C,  and  C„  from 

N  T 

such  data,  called  "bulk  methods".  In  particular,  a  method  was  suggested 
in  NEPRF  Tech.  Rept:  TR79-03,  based  largely  on  techniques  suggested  by 
Liu  et  al.  (1979).  Davidson  et  al.  (1981)  tested  this  method,  but.  only  over 
small  range  of  stability. 

We  discuss  in  this  report  an  attempt  to  test  the  NEPRF  bulk  method 
(with  slight  modifications)  on  a  more  complete  set  of  data,  obtained  by 
Donelan  on  a  mast  in  Lake  Ontario.  He  measured  fluctuations  of  velocity 


components,  moisture  and  temperature  as  well  as  profiles  of  these  variables 
during  a  wide  range  of  stability  conditions.  Surface  temperature  was  also 
available. 

The  purpose  of  the  study  was  three-fold. 

2 

1.  To  test  the  NEPRF  method  for  estimating  . 

2 

2.  To  evaluate  the  effect  of  moisture  on  C„  . 

N 

3.  To  test  bulk  methods  for  estimating  the  terms  contributing  to 

2  2 

C  and  C  ,  that  is,  for  estimating  vertical  flux  of  momentum, 

N  T 

heat  and  moisture.  These  quantities,  of  course,  have  uses  of 

2  2 

their  own,  independent  of  their  association  with  C_  and  C„  . 

T  N 

We  will  see  that  for  various  instrumental  reasons,  item  (2)  above 
was  not  resolved  satisfactorily;  but  the  other  two  items  were, 

2,  Observations 
2,1  General 

During  the  1972  International  Field  Year  for  the  Great  Lakes  (IFYGL) , 
experiments  were  conducted  to  measure  vertical  turbulent  fluxes  of  momentum, 
heat,  and  moisture  at  the  surface  of  Lake  Ontario.  Participants  in  these 
experiments  included  groups  from  the  Canadian  Center  for  Inland  Waters  (CCIW) , 
the  Atmospheric  Environment  Service,  Canada  (AES),  the  Bedford  Institute  of 
Oceanography,  and  many  other  groups  from  both  Canada  and  the  United  States. 

Mark  Donelan,  of  the  CCIW,  obtained  a  data  set  which  included 
fast-response  measurements  of  the  three  orthogonal-wind  components,  temperature 
and  humidity,  as  well  as  profile  measurements  of  the  mean  wind,  air  temperature 
humidity,  and  water  temperature.  These  observations  were  made  over  a  five- 
month  period  from  mid-May  to  mid-September,  thus  covering  a  wide  range  of 
stabilities.  It  is  this  data  set  which  is  used  to  test  a  modified  NEPRF 


bulk  model  for  the  estimation  of  structure  contrasts.  In  the  same  period, 
Canadian  meteorological  buoys,  scattered  over  Lake  Ontario  (Figure  la), 
recorded  air  and  water  temperature,  wind  speed,  direction,  and  humidity 
at  10  minute  intervals. 

2.2  Site  Description 

Donelan’s  observations  were  taken  at  a  location  2.6  km  off  the  south 
shore  of  Lake  Ontario  and  approximately  4  km  west  of  the  mouth  of  the 
Niagara  River  (Figure  la) .  The  sensor  masts  were  supported  by  guy  wires 
in  12  m  of  water.  A  barge,  located  approximately  170  m  southwest  of  the 
masts  (Figure  lb),  housed  the  recording  instruments  which  were  connected 
to  the  sensors  by  cables.  The  center  mast  in  Figure  lb  was  used  for  the 
turbulent  fluctuation  sensors  while  the  mast  50  m  to  the  east  held  the 
profile  sensors. 

This  site  offered  long  fetches  with  north  to  east  winds.  A  list  of 
the  wind  directions  and  their  associated  fetches  is  found  in  Table  1. 

Also,  the  nearness  to  shore  made  this  site  easily  accessible  by  small  boats. 

A  major  concern  about  this  location  is  its  proximity  to  the  mouth  of 
the  Niagara  River.  During  late  spring  through  early  autumn,  large 

horizontal  surface  temperature  gradients  are  created  where  the  water  of  the 
Niagara  River  empties  into  Lake  Ontario.  When  the  meandering  of  this 
"Niagara  plume"  brings  it  close  to  the  measurement  site,  inhomogeneities 
arise  in  the  surface  layer  due  to  the  development  of  an  internal  atmospheric 
boundary  layer.  Since  the  bulk  method  is  valid  only  under  homogeneous 
conditions,  cases  in  which  inhomogeneities  exist  will  give  unreliable 
results.  It  is  therefore,  important  to  identify  these  cases. 
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Fig.  1  (a)  Outline  of  Lake  Ontario  showing  the  positions  of  the 

meteorological  buoys  and  towers  during  IFYGL.  (From 
Donelan  et  al.  1974) .  (b)  Respective  positions  of  towers 
and  instrument  barge.  The  two  eastern  towers  were  used 
by  the  CCIW.  (From  Smith,  1974). 
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No  data  were  collected  during  periods  of  off-shore  winds  since  the 
air  Is  in  a  state  of  transition,  and  the  fetch  much  too  short  for  equi¬ 
librium  conditions. 

2 . 3  Instrumentation 

Table  2  briefly  describes  the  instrumentation  employed  for  the 
turbulent  fluctuation  measurements.  The  variables  were  all  sampled  at 
a  rate  of  20  times  per  second. 

Measurements  of  the  three  orthogonal  wind  components  were  made  at  a 
height  of  9.55  m  above  the  lake  surface  with  an  ultrasonic  anemomenter- 
thermometer.  The  instrument  was  a  model  PAT-313A-1  with  a  model  TR-31 
probe,  both  manufactured  by  Marine  Instrument  Company,  LTD.,  Tokyo,  Japan. 

The  measurement  consisted  of  a  pulse  travel  time  differencing  system  across 
a  probe  span  of  20  cm.  See  Mitsuta  (1966)  for  a  discussion  of  instrumental 
details.  The  instrument  had  a  resolution  of  +  2  cm  sec  1  and  frequency 
response  of  1000  Hz. 

Turbulent  fluctuations  of  the  temperature  were  measured  in  two  ways: 
by  the  same  anemometer-thermomenter  described  above;  and  a  platinum 
wire  sensor.  The  sonic  anemometer-thermometer,  at  a  height  of  9.55  m, 
employed  a  pulse  travel  time  summing  method  to  measure  the  temperature 
fluctuations  and  had  a  resolution  of  +  0.1°C  and  frequency  response  of 
100  Hz.  See  Mitsuta  (1966)  for  a  discussion  of  the  theory  of  the  techniques. 

The  cold  wire  instrumentation  consisted  of  a  model  1040  temperature 
and  switching  module  and  a  model  1210  standard  straight  platinum  wire 
probe,  both  manufactured  by  TSI  Incorporated,  St.  Paul,  Minnesota.  The 
wire  sensor  had  a  diameter  of  0.004  mm,  a  length  of  1.3  mm,  and  a  frequency 
response  of  5  k  Hz.  1  It  was  placed  at  4.95  m  and  later  moved  to  9.55  m  toward 
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Table  1.  Wind  directions  and  associated  fetch. 


direction 

S 

W 

NW 

N 

NE 


Table  2.  Sensors  used  for  turbulent  fluctuation  measureme 

sensor 

ultrasonic  anemometer 

ultrasonic  anemometer 

ultrasonic  thermometer; 
platinum  wire 

lyman-alpha  humidiometer 

dewpoint  hygrometer 


Table  3.  Sensors  used  for  profile  measurements 

sensor 

lightweight  cup  anemometer 
damped  wind  vane 
thermistor 

thermistor,  submerged 

standard  radiosonde 
carbon  hygristor 


the  end  of  the  experiment.  The  instrument  measured  the  ciectrica]  resistance 
of  the  wire,  which  is  directly  proportional  to  the  temperature  over  small 
temperature  ranges. 

The  turbulent  humidity  fluctuations  were  measured  by  a  lyman-alpha 
humidiometer  (LAH)  in  combination  with  a  thermoelectric  dewpoint  hygrometer. 

The  model  BL  lyman-alpha  humidiometer  used  was  manuf actured  by  Electromagnetic. 
Research  Corporation,  College,  Park,  Maryland,  and  had  a  frequency  response 
of  1500  Hz. 

The  LAH  Jias  a  sensor  tube  which  emits  a  beam  of  lyman-alpha  radiation 
(A  =  0.00122  cm)  through  a  lithium  fluoride  window,  across  a  measuring  path, 
through  a  second  window,  and  into  a  detector.  Due  to  the  lyman-alpha 
absorption  band  of  hydrogen,  ambient  water  vapor  in  the  measuring  path  decreases 
the  beam  intensity.  The  unabsorbed  portion  of  the  beam,  upon  entering  the 
detector,  ionizes  the  nitric  oxide  in  the  detector  and  generates  a  current 
according  to  the  relation 


where  1  is  the  generated  current,  c  is  the  constant,  x  is  the  path  length, 
and  a  is  the  absorption  coefficient  for  water  vapor  (380  cm  ^)  .  The  absorp¬ 
tion  coefficient  for  dry  air  is  the  order  of  0.34  cm  Therefore,  I  is 
essentially  a  measure  of  the  water  vapor  density,  q,  in  the  measuring  path. 
The  output  current  decreases  exponentially  with  an  increase  in  q,  but  a 
logarithmic  amplifier  was  used  to  produce  a  recorded  signal  varying  linerly 
with  q. 

Lithium  fluoride  was  used  in  the  sensor  as  window  material  due  to  its 
high  transmissivity  of  lyman-alpha  radiation.  A  disadvantage  of  this 


material  Is  that  It  Is  soluble  in  water,  so  that  sea  f; j*ray  or  condensation 
on  the  windows  will  etch  them  and,  therefore,  lower  their  transmissivity. 

For  this  reason,  a  dewpoint  hygrometer  was  used  In  conjunction  with  the 
LAH  for  constant  calibration. 

The  calibration  between  the  two  instruments  was  done  during  the  data 
analysis.  Also,  both  instruments  were  at  the  .same  height  for  all  the  runs 
analyzed.  Of  course,  the  response  of  the  hygrometer  is  too  slow  for  mea¬ 
surements  of  humidity  fluctuations  at  high  frequencies. 

The  dewpoint  hygrometer  system  used  was  a  model  137-C1,  Vapor  Mate  II 
with  a  model  137-C2  sensor,  both  manufactured  by  Cambridge  Systems,  Waltham, 
Massachusetts.  This  instrument  employed  a  cooled  mirror  which  was  warmed 
by  condensation  of  ambient  water  vapor.  A  precision  thermistor  measured 
the  temperature  change.  The  output  of  the  thermistor,  being  a  nonlinear 
function  of  temperature,  was  linearized  through  a  bridge  circuit. 

The  profiles  were  measured  on  a  tower  50  m  to  the  east  of  the  location 
of  the  measurements  of  the  turbulent  fluctuations.  The  wind  direction  was 
also  measured  at  a  height  of  15  m  above  the  lake  surface.  The  mean  wind 
speed,  temperature,  and  humidity  were  measured  at  fixed  levels  of  12.1  m, 

8.0  m,  5.2  m,  3. A  m,  and  2.1  m.  The  water  temperature  was  measured  a 
few  centimeters  below  the  water  surface  to  yie^d  a  ’’bucket  temperature".  A 
brief  description  of  the  sensors  used  for  the  profile  measurements  is 
found  in  Table  3.  Mean  values  of  the  measurements  at  the  start  and  end 
of  consecutive  10  minute  intervals  were  determined.  Donelan  et  al.  (1974) 
gives  a  detailed  description  of  the  sensors. 
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During  the  first  three  months  of  this  experiment,  profile  measure¬ 
ments  of  humidity  were  not  available  due  to  problems  with  the  carbon 
hygristor.  Instead,  mean  dewpoint  temperatures  from  the  dewpoint  hygro¬ 
meter  on  the  center  tower  at  9.55  m  were  used. 

There  were  also  eight  periods  in  which  turbulent  fluctuation  measure¬ 
ments  were  available  while  profile  measurements  were  not.  During  these 
periods,  bulk  methods  were  based  on  measurements  at  the  nearby  buoy  H, . 

2.4  Analysis 

The  turbulent  fluctuations  were  recorded  on  two  different  FM  tape 
recorders  in  analog  form  and  later  digitized  at  20  Hz  for  analysis  by 
computer.  Time  series  were  recorded  on  a  plotter  and  later  used  to 
identify  good  segments  of  data  to  be  used  as  individual  runs.  A  run  was 
usually  terminated  at  a  point  just  before  any  one  of  the  sensors  reached 
saturation  which  is  when  the  recorded  signal  from  a  sensor  goes  off-scale. 

Time  series  which  did  not  reach  saturation  were  usually  of  long  duration 

and  were  divided  into  two  or  more  separate  runs.  Table  4  shows  the  number, 
and  duration  of  each  run,  as  weii  as  special  characteristics  of  each  run. 

Once  the  runs  were  defined,  analysis  of  the  data  was  performed  on  a 
Cyber  171  computer  system  at  the  CCIW.  The  Cyber  171  utilized  a  fast- 
fourier  tranform  scheme  to  produce  spectra  cospectra,  and  quadrature 
spectra  of  the  three  orthogonal-wind  components,  sonic  temperature,  and 
humidity.  The  platinum  wire  temperature  fluctuations  could  not  be  corre¬ 
lated  with  any  of  the  sonic  or  humidity  fluctuations  since  they  were 
recorded  on  a  different  tape  recording  system.  Therefore,  only  one-dimensi on- 1 
temperature  spectra  were  obtained  from  the  wire  data,  and  the  instrument 
could  not  be  used  for  estimating  the  heat  flux. 
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The  coepectra  yielded  direct  measurements  of  the  fluxes  while  the 
one-dimensional  temperature  and  moisture  spectra  were  used  to  determine 
the  structure  constants. 

The  one-dimensional  horizontal-wind  spectra  were  also  used  to  determine 
the  momentum  flux  (stress)  through  the  dissipation  method.  The  different 
estimates  of  the  fluxes  and  structure  constants  will  be  compared  with  the 
bulk  estimates  in  Chaper  4. 

A  discussion  of  the  analysis  techniques  for  the  raw  profile  data  appears 
in  Donelan  et  al.  (1974).  Mean  values  of  wind  speed,  air  temperature,  and 
humidity  were  available  for  the  five  fixed  levels  mentioned  in  Section  2.3.2 
along  with  the  water  temperature.  As  was  mentioned  before,  humidity  profiles 
were  available  only  during  the  later  runs.  Since  the  NEPRF  bulk  model 
requires  the  mean  wind  at  10  m  for  estimating  roughness  lengths,  a  reference 
height  of  10  m  was  chosen.  Since  the  variables  were  usually  available  at 
8  and  12  m,  estimates  for  10  m  could  be  made  by  simple  interpolation. 

For  several  stable  runs  profile  information  was  not  available.  For  these 
cases,  a  reference  height  of  3.6  m  was  chosen  for  the  NEPRF  bulk  model  since 
this  corresponds  to  the  height  at  which  the  buoy  measurements  were  taken. 
However,  the  model  still  required  the  mean  wind  at  10  m  for  the  estimation  of 
roughness  lengths.  To  obtain  an  approximate  10  m  wind  speed,  the  4.0  m  wind 
was  multiplied  by  a  factor  of  1.2.  This  multiplication  factor  was  obtained 
from  the  average  ratio  of  10  m  to  4  m  winds  during  stable  periods  with  complete 
wind  information.  Although  this  factor  is  a  function  of  stability,  these  esti¬ 
mated  10  m  wind  should  give  satisfactory  estimates  of  the  roughness  lengths, 
which  do  not  vary  much  over  small  ranges  of  wind  speed. 

Tables  4  and  4a  list  the  runs  used  and  their  idiosyncrasies. 
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Kurt 

1972 

Date 

Date  Code 
Duration  (see 

(min)  Table  4a) 

Run 

1972 

Date 

Duration 

(min) 

Data  Code 
(see 

_ Table  4a) 

3 

May  22 

24 

A 

21-1 

27 

58 

B 

4 

22 

26 

A 

22 

June  16 

15 

B 

4-1 

22 

26 

A 

22-1 

16 

21 

B 

5 

24 

30 

A 

23 

16 

46 

A 

6 

24 

31 

A 

24 

19 

39 

B 

6-1 

24 

18 

A 

25 

21 

32 

A 

7 

24 

24 

A 

26 

22 

11 

A 

7-1 

24 

23 

A 

27 

24 

63 

A 

8 

25 

45 

A 

30 

Aug.  25 

90 

D 

9 

25 

13 

A 

31 

31 

82 

C 

10 

25 

53 

A 

32 

Oct.  5 

17 

C 

11 

25 

27 

A 

33 

12 

54 

E 

12 

25 

51 

A 

34 

13 

40 

E 

13 

25 

43 

A 

34-1 

13 

40 

E 

14 

25 

44 

A 

35 

13 

50 

E 

15 

26 

18 

A 

35-1 

13 

47 

E 

16 

26 

61 

A 

36 

13 

64 

E 

18 

26 

91 

A 

37 

13 

43 

E 

19 

27 

47 

B 

37-1 

13 

36 

- 

20 

27 

19 

B 

37-2 

13 

36 

E 

21 

27 

29 

B 

38 

13 

20 

E 

Table  4. 

Date  and 

duration  of 

each  run  used  for  analysis. 

and  data 

codes  defined  in  table  4a. 


Table  4a 


Data 

Code 


A 


B 


C 


D 


E 


Fast-Response  Sensor  Height 
4.95  m _  _ 9.55  m 


Cold  Wire 


Cold  Wire 


Ultrasonic 

LAH 


Ultrasonic 

LAH 


Ultrasonic 

LAH, 

Cold  Wire 


Ultrasonic, 

LAH, 

Cold  Wire 


LAH 

Cold  Wire 


Profile 

Information 

Available 


Yes 


No 


Yes 


No 


Yes 
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3.  Theory 

3.1  General  Characteristics  of 

Let  N  be  the  Index  of  refractions  for  the  wave  length  used  In  a  laser. 
Then  the  structure  function  for  this  index  is  defined  by  (x  is  a  position, 
variable) : 


D.,(r)  =  [N(x)  -  N(x  f  r)  V 
N 


Here  an  overbar  should  represent  an  ensemble  average.  In  practice  it  is 
taken  as  a  space  average. 


If  the  turbulence  is  homogeneous,  D„  depends  on  separation  r  only. 

N 

For  small  r  (r  <<  z,  the  height),  we  can  use  Kologorov  theory: 


n  ,  ,  -1/3  2/3 

DN(r)  =  ae  X ^ 


Here  "a"  is  a  constant,  about  3;  e  is  the  dissipation  of  turbulent  energy, 

2 

and  xN  is  the  rate  of  destruction  of  N*  /2  by  molecular  action  (where  N' 

denotes  the  deviation  from  an  average). 

-1/3 

The  quantity  ae  X^  is  defined  as  the  structure  "constant",  and 
2  2 

denoted  by  .  Given  ,  the  interaction  between  turbulence  and  laser 
beams  can  be  evaluated  by  the  "classical"  theory  derived  by  Tatarski  (1971). 
If  Monin-Obukhov  theory  is  adequate,  dimensional  analysis  gives: 


S2  =  N*2  z~2l\  (z/L) 


f 


where  g  is  a  universal  function.  is  defined  (in  analogy  with  T^)  by: 


N*  " 


w  N 


(2) 


Here  w  is  the  vertical  velocity,  and  =  -  w'T*  .  The  minus  sign  in  this 

u* 

definition  is  arbitrary,  and  added  so  that  T*  has  the  same  sign  as  Richardson 


number,  and  Monin-Obukhov  L,  defined  by  L  = 


cp0T  a  + 

kg  H  u  +  B  } 


-1 


where  k  is  the  vonKarman  constant  and  p  is  the  Bowen  ratio;  the  specific 
heat  at  constant  pressure;  g,  gravity  and  H,  the  vertical  heat  flux. 

Now,  N'  is  given  by  aT'  +  bq 1  where  a  and  b  are  known  parameters  which 
are  constant  for  practical  purposes,  q  is  the  specific  humidity.  Hence, 
equ.  (2)  becomes 


N .  *  a 
*  u 


+  b  =  -  aT*  -  bQ, 

*  * 


(3) 


where  is  defined  as  -  w'q'/u^. 
Hence , 


N*  =  +  a  TJ1  +  Y  T  ] 
1 * 


where  y  m  b/a. 

Q* 


But  —  -  C _/l& 


where  6  is  the  Bowen  ratio  and  Z  the  latent  heat. 


(4) 


Hence , 


2 

is  defined  through  the  structure  function  and  Kolmogorov  theory 

2 

in  analogy  with  the  definition  of  .  Therefore,  we  can  put  (see  Wyngaard, 
et  al. ,  1971) . 


f (z/L) 


(6) 


If  f(z/L)  is  the  same  universal  function  as  g(z/L)  above,  equations 
(1)  and  (5)  combine  to  form 


2  2 

a^CT  [1  + 


L1  ) 


(7) 


Y  C 

The  ratio  _ £. 


is  0.03  for  lasers. 


I 

Hence,  finally 


[1  + 


0.03.2 
B  J 


(8) 


2 

If  this  equation  were  correct,  only  estimates  of  CT  and  Bowen  ratio 

2 

would  be  needed  to  obtain  C„  .  However,  equ.  8  is  controversial.  Mcllveen 

N 

(1981)  derived  instead  (see  also  Wesely,  1976) 


u  + 


.06 

6 


r  + 


(M) 
V  6  ’ 


(9) 


Where  r  is  the  correlation  coefficient  between  moisture  differences  and 
temperature  differences  at  points  with  small  separations.  If  r  *  1,  or  -1 , 
equations  8  and  9  coincide,  since  r  will  have  the  same  sign  as  6. 

We  had  hoped  to  test  the  correlations  between  temperature  fluctuations 


and  moisture  fluctuations  in  this  project.  But  only  sonic  temperatures 


were  measured  at  the  same  times  and  nearly  the  same  place  as  moisture. 

And,  as  we  will  see,  the  sonic  temperature  spectra  were  extremely 

2 

noi6y  and  unreliable.  We  expect  r  to  be  close  to  1,  so  that  the  difference 
between  equations  8  and  9  is  not  important. 

2 

In  most  cases,  the  influence  of  moisture  on  C„  is  expected  to  be  small. 

N 

In  only  2  of  the  41  cases  for  which  bulk  Bowen  ratios  were  computed, 

the  absolute  value  of  8  was  less  than  0.1.  Also,  in  2  cases,  { 3 |  exceeded 
1.  Hence,  usually,  the  difference  between  equations  8  and  9  is  not  of  critical 
importance,  over  Lake  Ontario.  This  does  not  mean  that  the  effect  of  moisture 

is  always  unimportant  over  Lake  Ontario:  in  the  two  cases  in  which  B  was  less 

2 

than  0.1,  the  moisture  effect  should  not  be  neglected  in  the  evaluation  of  C„  . 

N 

2 

But,  in  these  cases,  CT  was  quite  small. 

2 

Here,  we  will  concentrate  here  on  CT  .  We  will  begin  with  equ  (6),  written 
in  the  form 


«  ( 


u* 


2 

) 


g(z/L)z  2^3 


(10) 


As  we  have  seen,  the  length  L  depends  only  on  w’T' ,  u^  =  ~u  V  ,  and  8, 
which  also  depends  on  w'q'  and  w'T’. 

Hence,  we  will  first  consider  bulk  methods  to  determine  w’T',  w'u'. 
and  w’q’,  from  ship-board  data. 

3.2  Bulk  methods  for  stress 

For  aerodynamic  flow,  u’w'  is  usually  described  by: 


-  u'w'  = 


.2  2 

k  u 

z 


dnf-  -  Hh) 

Zo  L 


(11) 


where  u  is  the  wind  speed  at  height  z,  z  is  the  roughness  length,  and  $ 
z  o 


is  a  presumably  universal  function  which  is  quite  well  known  and  has  been 
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tabulated.  It  is  negative  in  stable  air  and  positive  in  unstable  air. 

2  2 

The  factor  k  /(In  z/z  -  ’■Kt-))'  is  often  called  the  "drag  coefficient",  T  . 

o  L  m 

There  is  some  disagreement  as  to  the  choice  of  z  .  It  is  usually 

o 

between  0.01  and  0.5  cm  and  increases  with  increasing  wind  speed.  Dif¬ 
ferent  investigators  have  derived  somewhat  different  relationships 

between  z  and  wind  speed.  In  this  work,  we  have  made  use  of  three 
o 

different  relations,  suggested,  respectively tby  Liu  et  al.  (1979), 

Smith  and  Banke  (1975)  and  Garratt  (1977).  The  differences  between  these 

three  formulations  are  not  very  great.  More  serious  is  the  problem  that 
the  stress,  -  P  uTrcr  ^  should  not  only  depend  on  the  wind  speed, 

roughness  and  stability,  but  also  on  the  difference  between  wind  speed 
and  wave  speed — which  depends  on  the  "age"  of  the  wave;  in  well-developed 
waves  with  large  fetch,  wave  speed  and  wave  speed  can  be  very  similar, 
reducing  the  stress  below  that  expected  from  equation  (11). 

In  order  to  evaluate  the  stress  by  equ  (11),  we  employ  a  process  of 

successive  approximations,  which  converges  rapidly  in  practice.  We  first 

estimate  u’w'  for  neutral  conditions,  and  use  the  corresponding  technique 

to  infer  w'T'  from  a; r-sea  temperature  difference  and  temperature  roughness 

length  for  =  0  (see  next  section).  Given  w'T'  and  u'w',  we  compute 

z  - 

L  and  hence  and  recompute  u'w'  and  L.  This  procedure  usually  has  to 

la 

be  repeated  only  twice  to  yield  stable  estimates  of  u'w'. 

3.3  Bulk  method  for  heat  flux 

For  heat  flux,  Monin-Obukhov  theory  suggests: 
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where  0  Is  the  (potential)  temperature  at  height  z  and  G  is  the  surface 
A  s 

temperature.  Again,  the  function  i^(j-)  is  well  known.  One  difficulty  with 
the  use  of  (12)  is  that  0^  is  the  "skin"  temperature,  whereas  "bucket" 
temperatures  are  usually  measured.  The  bucket  temperature  can  be  signifi¬ 
cantly  warmer  than  the  skin  temperature 

In  order  to  correct  for  the  possible  effect  of  tiiis  difference,  Liu 
et  al.  (1979)  have  constructed  a  nomogram  which  permits  estimates  of  differ¬ 
ences  between  skin  and  bucket  temperatures  from  wind  speed  and  differences 
between  air  and  bucket  temperatures.  We  use  the  results  of  this  nomogram  in 
the  numerical  form  suggested  by  Burk  et  al.  (1979). 

Another  problem  is  the  interpretation  of  z  ,  This  is  not  really  a  func¬ 
tion  of  ocean  roughness,  but  a  measure  of  the  molecular  transfer  of  heat 
through  the  surface.  It  is  of  order  . /u^  where  is  the  molecular  heat  con¬ 
ductivity.  We  use  here  the  estimate  of  z  based  on  Liu  et  al.'s  relations 
under  neutral  conditions  as  functions  of  wind  speed. 

We  then  proceed  to  soLve  equ  12  first  for  ..  =  ■).  This  is  combined,  as  be¬ 
ll 

fore,  with  -u’w’  computed  for  z/L  =  0,  to  form  1.  Hence  the  computations  of 
w’T’  is  repeated  several  times  leading  to  presumably  improved  estimates  of  L. 

3.4  Evaporation 

Although  w’q’  is  not  usually  important  for  ,  bulk  estimates  ot 

evaporation  are  important  for  other  reasons.  Further  the  effect  of  w’q’  on 

C  2  cannot  be  neglected  over  tropical  oceans.  In  analogy  with  equations 
N 

(11)  and  (12),  we  then  have: 


w'q ' 


k“(q  -q  )u 
s  z  z 


[In 


oq 


Vl)]  [ln  r  -  *<t? 


(13) 
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Here,  qg ,  the  surface  specific  humidity,  is  evaluated  from  the  surface 

skin  temperature;  4^  is  equated  with  4^  and  z  is  also  obtained  from 

Liu  et  al.  (1979).  Then,  equ.  13  is  solved  by  successive  approximations, 

initially  with  L  =  ®,  and  then  with  progressively  more  realistic  L-values. 

This  whole  procedure  is  analogous  to  that  suggested  by  Burk  et  al. 

(1979)  except  that  z  ,  z  _  and  z  are  assumed  to  depend  on  wind,  and  not 

o  oT  oq 

on  stability.  Thus,  these  parameters  are  not  changed  in  successive 
iterations . 


3. 5  Other  techniques  for  estimating  fluxes  and  structure  constants 

The  quantities  -u'w'  could  be  estimated  directly  from  observations 
of  fluctuations  of  u  and  w.  However,  the  results  are  extremely  sensitive 
to  errors  in  tilt  of  the  mast. 

Another  technique,  depending  on  accurate  high-frequencies  recording 
(frequency  >  1  hz)  is  based  on  Kolmogorov  theory: 


k1S(k1)  =  c 


2/3  .f. 
£  <u> 


-2/3 


(14) 


where  S  is  the  spectral  density  at  longitudinal  wave  number  k^,  c  is  a 
universal  constant,  f,  frequency,  and  c,  as  before,  energy  dissipation 
into  heat. 

For  Monin-Obukhov  scaling,  e  is  described  by: 


e 


(15) 
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We  use  here  <t>^  =  1  -  ^  in  unstable  air,  and  =  [1  +  2.5  (^O'^]^2  in 

stable  air.  Hence,  if  we  eliminate  e  between  equ.  (15)  and  (14)  we 

2  - 

derive  an  equation  for  u^  =  -  u'w  . 

w'q'  and  w'T'  were  measured  from  their  definitions  only,  although 
there  exist  methods  to  infer  them  from  high-frequency  fluctuations  of 
q  and  T,  respectively,  once  c  or  u  are  given. 

2 

CT  was  estimated  from  high-frequency  temperature  spectra  (k^  in 
radians  per  unit  length) : 


ST(k1)  =  1/4  CT2  kx  "5/3 


2mf 


with  Taylor's  hypothesis  that  =  — —  (f  in  hz) . 


An  attempt  was  also  made  to  estimate  for  sound  from  the  analogous 

equations  for  index  of  refraction  for  sound: 


W  "  i  S2  \ 


-5/3 


making  use  of  the  spectra  of  the  sonic  thermometers,  in  order  to  check 
equ.  (9). 
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■■  n?* 


4. 


Results 


4 . 1  Heat  f lux 

Figure  2  compares  directly  measured  w'T'  with  w'T'  estimated  from  the 
ulk  method  (equ.  12).  Since  onlv  sonic  temperatures  could  be  used,  the 
unstable  cases  of  October  are  not  included.  The  sonic  temperatures  had  been 
corrected  for  moisture  and  wind  effects.  The  general  agreement  is  surpris¬ 
ingly  good.  There  are  a  few  cases  with  generally  small  heat  flux,  when  the 
directly  measured  Reynolds  flux  was  negative,  but  the  bulk  heat  flux  was  posi¬ 
tive.  For  some  of  these  cases,  temperature  profiles  were  available.  These 
showed  that  the  air  was  slightly  stable  around  10  m  where  w'T’  was  measured, 
in  agreement  with  the  negative  sign  of  w'T'.  But  the  surface  temperature  was 
larger  than  the  temperature  at  10  m  leading  to  upward  heat  flux  by  the  bulk 
method.  This  kind  of  profile  is  not  homogeneous  but  suggests  that  the  lowest 
layers  had  been  warmed,  perhaps  when  the  air  drifted  over  the  near-by  warm 
Niagara  plume. 

If  the  bulk  method  value  of  flux  had  been  computed  from  a  surface  tempera¬ 
ture,  extrapolated  downward  from  the  profile,  the  bulk  fluxes  would  also  have 
been  negative  and  in  good  agreement  with  measured  w'T’.  For  large  downward 
flux,  there  appears  a  tendency  for  the  bulk  flux  to  be  larger  in  absolute 

value  than  w'T'.  On  the  average  there  is  a  slight  bias  in  the  sense  that  bulk 

_2 

heat  fluxes  exceed  Reynolds  heat  fluxes  by  about  7  watts  m  “.  The  reason  lor 
this  bias  is  not  known.  Both  the  corrections  for  wind  and  moisture  fluxtuations 
are  in  the  wrong  direction,  so  that  tiie  agreement  could  have  been  improved 
slightly  if  these  corrections  had  not  been  made. 

4 . .  Stress 

Figure  3  shows  estimates  of  t./r  =  -  u'w'.  Bulk  results  are  shown  along 
the  abscissa,  direct  Reynolds  stress  measurements  and  c  results  along  the 
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u*  Bulk  (ro2/sec2) 


Fig.  3  Comparison  of  momentum  fluxes  from  eddy  correlation  method 
(dots)  and  dissipation  method  (circles)  with  bulk  metnou. 

! 
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ordinate.  Again,  Reynolds  stresses  and  dissipation  stresses  had  to  be  based 
on  sonic  anemometers,  so  that  the  unstable  cases  in  October  are  not  included 
in  this  figure.  For  a  fixed  bulk  estimate,  the  observed  stresses  by  dissipation 
and  direct  methods  are  not  usually  far  apart.  We  therefore  suggest  that 

the  best  estimates  of  "true"  stresses  are  close  to  these  measurements. 

2  2-2 

For  u^  <  0.04  m  sec  ,  the  agreement  between  bulk  methods  and  the 
others  is  reasonably  good.  Even  then,  the  bulk  estimates  (derived  from  the 
work  of  Liu,  et  al.  (1979)  tend  to  be  relatively  large.  But  in  cases  of 
large  stresses  (strong  winds)  the  bulk  methods  seem  to  overestimate  stresses 
by  about  a  factor  of  2  (or  u*  by  ,  2) . 

The  reason  is  not  clear.  Possibly,  in  a  well-developed  wave  pattern 
with  large  fetch,  wave  speeds  in  strong  winds  are  not  very  different  from 
wind  speed,  thus  producing  a  small  stress  for  a  given  wind. 

To  test  this  hypothesis,  we  plotted  the  differences  between  bulk  stress' 
and  Reynolds  stress  as  function  of  wind  direction  and  speed;  the  fetch  was 
a  strong  function  of  wind  direction,  with  longest  fetches  with  North  winds 
and  shortest  for  South  winds  (see  Table  1) . 

Figure  4  shows  the  results.  Indeed,  largest  differences  between  bulk 
and  Reynolds  stresses  occur  with  strong  winds  at  10  m  with  large  fetch,  sug¬ 
gesting  that  small  drag  coefficients  are  appropriate  in  these  cases  because  of 
the  small  differences  between  wind  speed  and  wave  speed.  Note  that  these  strong 
winds  were  directly  measured,  not  extrapolated  vertically  from  4  m. 

Since,  apparently,  bulk  estimates  overestimate  u^,  they  should  also  seri- 

3 

ously  overestimate  L,  which  is  proportional  to  u^  . 

4. 3  Moisture  flux 

Figure  5  compares  bulk  and  Reynolds  moisture  flux.  The  result  is  quite 
satisfactory  with  no  strong  bias  suggested.  There  is  a  slight  tendency,  for 
the  bulk  methods  to  overestimate  the  fluxes  when  w’q’  <  1  -  ,  .  For 

car*  Ir  o 


24 


E  gu|k  (g/m2  sec) 


Comparison  of  moisture  fluxes  from  eddy-correlation  method 
and  bulk  method.  Stable  periods  are  represented  hy  dots; 
unstable  by  squares;  neutral  by  triangles;  and  cases  when 
the  observed  and  bulk  sensible  heat  fluxes  have  opposite 
signs  by  open  circles. 


larger  fluxes,  no  bias  is  obvious.  However,  there  is  considerable  random 
scatter.  Note,  that,  again  the  unstable  cases  of  October  are  excluded. 


4.4  CT2  and  C^2 

—————  2  2 

Figure  6  compares  bulk  estimates  of  with  estimated  from  tem- 

2  -3  2 

perature  spectra.  The  agreement  is  especially  good  when  >  2  x  10  °C 

-2/3  2 

z  ,  both  for  stable  and  unstable  air.  For  very  small  ,  the  spectral 

estimates  are  significantly  higher  than  the  bulk  estimates.  We  believe  this 
is  due  to  noise  in  the  spectral  data;  since  spectra  are  positive  definite, 
any  "noise"  will  produce  net  positive  errors.  Also,  some  of  the  worst  dis¬ 
crepancies  occur  for  cases  where  the  bulk  heat  flux  and  measured  heat  flux 
had  opposite  sign,  so  that  the  basic  conditions  for  bulk  methods  were  not 
satisfied. 

2  2 

The  good  agreement  between  "observed"  CT  and  bulk  may  appear  sur¬ 
prising  in  view  of  the  poor  performance  of  bulk  methods  for  estimating  stress. 
The  reasons  may  be  these: 

w  *  T  * 

First,  the  determination  of  the  quantity  T*  =  - - —  by  bulk 

methods  is  independent  of  the  drag  coefficient  of  momentum.  This  is  because 
the  temperature  profile  is  given  by: 


_z 

L 


The  bulk  method  determines  T^  from  this  equation.  The  drag  coefficients 

2 

comes  in  only  weakly  through  <!>,  —  ,  since  L  depends  on  stress.  But  because 

n  (LJ 

the  incorrect  drag  coefficients  occur  with  strong  winds,  ^  is  small  compared 

to  In  z/z  _  in  these  cases. 
oT 

2 

The  drag  coefficient  enters  the  estimation  of  CT  also  through  the  func¬ 
tion  g(z/L)  in  equ.  10.  But  again,  the  doubtful  drag  coefficients  occur  f  m 
strong  winds,  where  g(z/L)  ^  g(0)  which  does  not  depend  on  stress. 
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Fig.  6  Comparison  of  temperature  structure  constants  from  one¬ 
dimensional  temperature  spectra  and  the  bulk  model.  Dots, 
stable;  squares,  unstable;  triangles,  neutral;  circles, 
bulk  and  Reynolds  fluxes  have  opposite  sign. 
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There  probably  is  little  advantage  in  applying  moisture  corrections 

2 

as  suggested  by  equations  (8)  or  (9)  in  order  to  calculate  C  ,  for  obser- 

N 

vations  such  as  those  available  here.  The  reason  is  that  the  corrections 
depends  on  the  Bowen  ratio.  But  here,  any  Bowen  ratio  estimated  by  bulk 
methods  may  not  even  agree  in  sign  with  the  actual  Bowen  ratio,  when  con¬ 
ditions  are  not  homogeneous.  E.g.  sometimes,  when  Tg  >  T,  the  bulk  heat 

flux  is  upward,  yet  the  local  heat  flux  at  10  m  is  down.  Moisture  flux  is 

upward,  in  all  cases,  so  that  the  bulk  Bowen  ratio  has  the  wrong  sign. 

4.5  Comparison  of  sonic  and  wire  temperature  spectra 

Sonic  thermometers  actually  measure  spectra  of  the  index  of  refraction 

for  sound.  This  is  proportional  to  the  speed  of  sound  given  by  20.1 '/i'O-  +  .51q). 

A  fluctuation  of  T  ,  T  ’ ,  is  then  approximately  given  by 
s  s 

T  *  -  Tr  +  .51  T  q' 
s 

If  the  correlation  between  q’  and  T'  is  high,  in  absolute  value,  spectral 
densities  are  releted  by 

_  2 

ST  =  ST  +  (.51T)  Sq  +  1.01  T  SqST 
s 

In  most  cases,  the  terms  involving  fluctuations  of  q  were  small,  so 

that  ST  should  be  nearly  equal  to  ST.  Fig.  7  compares  typical  spectra 
8 

of  T  and  T  for  a  case  of  small  moisture  variation.  Clearly,  the  agreement 
s 

is  poor,  especially  in  the  inertial  range,  where  the  sonic  spectra  remain 
high  when  they  should  decay  rapidly  with  frequency.  This  difference  is 
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fig.  7  Comparison  of  one-dimensional  temperature  spectra  for 
Run  21-1  from  cold  wire  (dots)  and  sonic  thermometer 
(circles),  w'q'  for  this  run  was  small. 
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most  likely  caused  by  "noise"  in  the  sonic  data,  and  has  the  consequence 
2  2 

that  ,  and  therefore,  for  sound  could  not  be  determined  from  the  sonic 
spectra . 

The  hypotheses  that  high-frequency  variations  of  sonic  temperature  are 

noisy  is  confirmed  by  the  chaotic  behavior  of  high  frequencies  of  the 

cospectra  between  sonic  temperature  and  moisture.  Therefore,  we  could  not 
2  2 

evaluate  (for  sound)  independently  from  ,  to  test  equations  (6)  or 
(9). 

Another  question  associated  with  these  equations  is  the  correlation 
between  high-frequency  variations  of  moisture  and  temperature.  If  this 
correlation  is  +  1  or  -  1,  these  equations  become  identical. 

Observations  indeed  show  numerically  large  correlations  between  temperture 
and  moisture  which,  however,  decrease  toward  high  frequencies.  This  may 
be  because  moisture  and  temperature  are  not  measured  at  the  same  location. 

It  seems  intuitively  reasonable  that  high  frequency  temperature  and  moisture 
fluctuations  at  a  fixed  location  should  be  well  correlated  since  they  are 
presumably  produced  by  the  same  high-frequency  vertical  velocities. 


4.6  C 

q 

2  2 

Figure  8  compares  bulk  estimates  of  C  with  C  determined  from 

q  q 

p 

moisture  spectra.  In  stable  air,  the  bulk  methods  underestimate  the  C 

q 

by  an  average  of  2,  and  in  unstable  air,  bulk  methods  tend  to  overestimate 
2 

the  Cq  by  an  average  factor  of  1.7,  at  least  for  cases  in  which  bulk 
heat  flux  and  local  heat  flux  at  10  m  have  the  same  sign.  Also,  in  unstable 
air,  the  scatter  is  quite  large.  But  on  the  average,  the  agreement  is 
reasonably  good. 
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II. 


We  can  write  C 

q 


in  analogy  with 


9 

C 

q 


-2/3 

7.  g(z/l.) 


In  order  to  try  to  understand  the  different  effect  of  stability  on 
9  9 

bulk  Cq“  and  spectral  ,  we  have  further  analyzed  the  behavior  of  g(z/L). 

We  had  previously  suggested  that  g(z/L)  is  the  same  funcation  as 

f(z/L)  in  equ.  6.  We  have  now  recomputed  g(z/L)  from  the  spectral  obser-  i 

2  ! 

vations  of  C  ,  and  show  it  in  Fig.  9  as  function  of  z/L.  Apparentlv,  J 

q  t  : 

it  differs  significantly  from  Wvngaard  et  al.'s  (1971)  f(z/L);  and  it 

differs  in  the  opposite  sense  as  f(z/L)  computed  by  Fairall  et  al.  (1980). 

3 

This  difference  may  be  due  to  the  overestimate  of  u^  by  the  bulk  method 
noted  earlier,  resulting  in  an  overestimate  in  the  absolute  value  of  I..  This 
means  that  the  bulk  method  underestimates  jz/hf.  If  we  increased  |z/l|  by  a 
factor  of  3  (as  suggested  by  fig.  3)  the  agreement  between  Wvngaard's 
equation  and  the  curve  inferred  from  the  present  data  would  be  quite  good. 

However,  we  must  remember  that  the  stress  comparison  excluded  most  of  the 
unstable  cases. 


1 

i 


33 


Fig.  9 


g(— )  according  to  Wyngaard 
L  2  2 

observed  C  “  and  bulk  Q.“; 
q 


(solid  line) 
second  solid 


dots  from 
line  drawn 


3  4 


to  dots. 


5.  Summary 


Observations  of  means  and  fluctuations  of  wind  components  and  tempera¬ 
ture  on  masts  over  Lake  Ontario  have  been  used  to  evaluate  bulk  methods  for 
2 

estimating  0^“,  and  fluxes  of  momentum  heat  and  moisture.  The  observations 

•xtended  from  generally  stable  conditions  in  May  to  unstable  conditions  in 

Oc  tober . 

2 

C  “  was  directly  determined  from  temperature  spectra  in  the  inertial 
range;  the  fluxes  were  estimated  by  the  product-momentum  method.  In  addition, 
stresses  could  be  determined  from  the  spectra  of  the  longitudinal  wind  com¬ 
ponents  in  the  inertial  range. 

Stresses  determined  by  the  two  methods  based  on  wind  fluctuations  agreed 
quite  well  with  each  other,  but  stresses  estimated  from  mean  winds  and 
drag  coefficients  were  significantly  larger  for  strong  winds  and  large  fetches. 
In  other  words,  the  drag  coefficients  assumed  from  previous  work  were  much 
too  large. 

Observed  and  bulk  heat  fluxes  agreed  quite  well.  There  were,  however, 
a  few  cases  of  weak  heat  flux  when  the  signs  derived  from  the  two  methods  were 
inconsistent.  In  these  cases,  the  local  lapse  rate  at  10  m  was  stable  (hence 
w'T'  <  0);  but  the  lake  surface  temperatures  were  larger  than  the  temperatures 
at  10  m  (hence  upward  bulk  flux) .  Here  the  vertical  structure  was  nonhomo- 
geneous,  probably  due  to  the  relatively  warm  Niagara  plume;  in  these  cases, 
the  conditions  for  use  of  the  bulk  methods  were  violated. 

Bulk  and  direct  estimates  of  moisture  flux  agreed  quite  well,  however, 
with  considerable  random  scatter. 

2 

Directly  measured  and  bulk  methods  for  estimation  of  CT  agreed  sur- 

2 

prisingly  well,  except  in  a  few  cases  of  small  C^,  where  the  signs  of  the  huix 
heat  fluxes  were  incorrect. 
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The  agreement  was  good  even  for  strong  winds  and  large  fetch,  when 

the  drag  coefficients  used  seriously  overestimated  the  stresses. 

It  is  shown  that,  in  strong  winds,  drag  coefficients  (for  momentum) 

2 

have  little  effect  on  . 

2 

C„  could  not  be  measured  directly;  but  the  measurement  of  the  fluxes 
N 

over  Lake  Ontario  and  theoretical  argument  suggest  that  moisture  only  rarely 

2  2  2 
affected  ,  so  that  was  essentially  proportional  to  C^,  at  this  site. 
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